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Abstract 

We develop the canonical ADM approach to 2+1 dimensional gravity in presence 
of point particles. The instantaneous York gauge can be applied for open universes or 
universes with the topology of the sphere. The sequence of canonical ADM equations 
is solved in terms of the conformal factor. A simple derivation is given for the solution 
of the two body problem. A geometrical characterization is given for the apparent 
singularities occurring in the N-body problem and it is shown how the Garnier hamil- 
tonian system arises in the ADM treatment by considering the time development of 
the conformal factor at the locations where the extrinsic curvature tensor vanishes. 

1 Introduction 

After the seminal paper by Deser Jackiw and 't Hooft the problem of supplying a solution 
to 2 + 1 dimensional gravity has been approached from different viewpoints and with different 
techniques. Moncrief [0 and Hosoya and Nakao gave the hamiltonian treatment in absence 
of particles. In [0, |^ the complete reduction of the hamiltonian to the physical parameters i.e. 
the moduli of the time slices was obtained in the York (minimal surface) gauge; the problem 
can be dealt with explicitly for the torus while for higher genus, even if well defined, is far 
from trivial. The solution of and [Q] were exploited in 0] and for the quantization of 
the theory. 

In presence of particles progress was made in the papers by Bellini, Ciafaloni and Val- 
tancoli and by Welling |^ in the first order formalism by going over to the instantaneous 
gauge; in these papers it was shown that the problem is equivalent to the solution of a 
Riemann- Hilbert problem. It can be solved explicitly for the case of two particles, in terms 
of hypergeometric functions. For three or more particles one encounters a feature known 
in the mathematical literature as apparent singularities. Such apparent singularities evolve 
according to a hamiltonian system of equations known as the Garnier system 0] which is 
derived in the mathematical literature from the isomonodromicity condition. 

A different approach was put forward by 't Hooft by describing the evolving Cauchy 
surfaces in terms of polygonal tiles which join along segments where the extrinsic curvature 
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is singular. The dynamics of the system is codified in transition rules which intervene when 
e.g. the length of a side of a polygon goes to zero or when a particle collides with the side 
of a polygon. A similar but different approach was given by Waelbroek |TD[. Quantization 
schemes were given in [0, 0, |l^ . 

In this paper we shall consider 2 + 1 dimensional gravity in presence of massive particles 
by exploiting the hamiltonian formulation. Most techniques apply to massless particles as 
well, but here we shall not be concerned with this case. Thus our setting is a second order 
formalism and the basic equations are those derived systematically from the variation of 
the ADM action. The reason for such a development is to give a treatment which resides 
completely within the canonical framework. We shall see that all results obtained in p 
[0 come out in simple fashion from such an approach; moreover we shall prove that the 
Garnier hamiltonian system is a direct outcome of the canonical ADM formalism. Again 
the exploitation of the instantaneous York gauge plays a major role; the technical advantage 
of such a gauge is to reduce an equation of sinh-Gordon to one of Liouville type to which 
powerful methods of complex analysis apply. Unfortunately the applicability of such a gauge 
is restricted to universes of spherical topology or to open ones. 

We shall see that such an approach provides an elementary way to solve the two body 
problem: the exact solution of the motion can be derived even before the explicit computation 
of the metric, which combined with the qualitative knowledge to the asymptotic metric gives 
a complete description of the scattering. The exact metric is obtained by solving a Liouville 
equation, leading of course to the same results obtained in |0 and 

One of the aims of presenting a canonical solution of the problem is to provide a framework 
for the quantization of the theory but here we shall be concerned solely with the classical 
theory. 

The plan of the paper is the following: After recalling in sect. 2 the ADM canonical for- 
mulation of gravity coupled with point particles we proceed in sect. 3 to the solution of the 
canonical equations. These will be solved in the following order: First one computes the 
momenta canonically conjugate to the space metric; they turn out to be rational analytic or 
antianalytic functions of the complex coordinates of the plane. The knowledge of such canon- 
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ical momenta allows one to write down the partial differential equation for the conformal 
factor or better for a more fundamental quantity which we shall call the reduced conformal 
factor. This is a Liouville equation whose sources are located at the particle position; in 
addition when the number of particles H is greater than 2 new sources appear. These are 
located at the points where the extrinsic curvature tensor vanishes and are identified as the 
position of the apparent singularities of the fuchsian differential equation which underlies 
the expression of the conformal factor. In the ADM equation for the conformal factor a the 
sources appear as function of the conformal factor itself. On the other hand it is shown in 
the text that due to the vanishing of the inverse g'^^ of the spacial metric on the particles, the 
sources simplify and arc just given by the particle masses. The Liouville equation admits a 
whole family of solutions parametrized by a number M whose values runs from zero to Aix 
and it is related to the asymptotic behavior of the conformal factor; it is nothing else than 
the total energy of the system in adimensional units. A byproduct of the ADM equation for 
the time derivative of a is that such parameter M is constant i.e. the energy is conserved. 

The Lagrange multiplier N which appears in the ADM metric in the combination —N'^dt^ 
satisfies a differential equation, which again due to the vanishing of g^^ on the particles 
reduces to a homogeneous linear differential equation. This is shown to be simply solved 
by the derivative of the conformal factor with respect to the parameter M. An elementary 
theorem proves that such a solution is unique. The final ingredient of the ADM metric i.e. 
the shift functions N'^ are expressed in terms of the spacial derivative of N. 

In sect. 4 we write the Hamilton equation for the time variation of the particle positions 
and momenta. These are explicited for the general A/'-body problem in terms of the coefficient 
of the fuchsian equation which underly the conformal factor. In the two body case, which 
is dealt with in sect. 5, they give rise to an elementary system of two differential equations 
for the particle positions and momenta, which solved provides the particle trajectories. The 
exact metric can be given in term of the classical solution of the Liouville equation. In sect. 6 
we give the general form of the asymptotic metric in the York instantaneous gauge. 

In sect. 7 we give the treatment of the problem for jV" > 2. It is shown how the fuchsian 
nature of the differential equation is related to the conservation of momenta and to a gen- 
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eralized conservation law which appears related to the transformation of the action under 
space dilatations. 

As mentioned above, with more than two particles there are apparent singularities and 
one has to provide evolution equations for their positions and residues. We show that 
all this information in contained in the ADM equation for the the time derivative of the 
conformal factor. In fact by exploiting Schwarz's relation between the coefficient of the 
fuchsian differential equation and the reduced conformal factor one obtains the Garnier 
equation by equating the residues at the polar singularities at the position of the apparent 
singularities. 

2 The action 

The action for the gravitational field in any dimension, including the boundary terms is 



given by [0 



„3 r 3 j-Bi 

oTTGat Jb SttGn J Bo 

where g^u, are the -D + 1 dimensional metric and curvature, gij the D-dimensional metric 
of the constant time slices S^; Sq and Ei the initial and final time slices; Kij the second 
fundamental form of the time slices and 6jj the second fundamental form of the lateral 
boundary B whose volume form is v^— 7; is the future pointing unit normal to the time 
slices and the outward pointing unit normal to B; Bt = Cl B, ^Ja is the volume form 
induced on Bt and sinhi] = n^u^. In presence of particles we must add to the gravitational 
action the matter term 



j dt iRnC^ -g^^iqn)qnqn (2) 
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where we have chosen the gauge g° = t. At the classical level it is convenient to multiply 
the action by ^^^r^; the matter action e.g. becomes 

Sm = - dt^nin^J -g^y{qn)qnq';, (3) 

with nin = i^l^^; it is a real number < m„ < An. In the ADM notation where the 
metric is written as ||T^ 

ds'^ = -N'^dt^ + Qij {dx' + N'dt) {dx^ + N^dt) (4) 

the volume term can be rewritten as [ITH 



Sh= I d^xdt N^[R + KijK'i - + 2{n^K - a^").^] , (5) 

J M 

where R is the D dimensional curvature related to gij and is the acceleration vector 
= n^n^.j^. Going over to the canonical variables 

Fni = -wr- = m„ = = m„ — = (d) 

the matter action can be rewritten as 

Sm= jdtY, [Pmcfn + N'^q^)?^, - N {q„) ^^^^^9^0^+^) ■ (7) 
On the other hand the gravitational action expressed in terms of the canonical variables 



becomes |Tj, |I6|, [17 



Scrav = I dtd^'x [n'^g,, - N'H, - NH] + 
Jm 

+ 2 jdt j d(^-i)a; v^iV (^i^B, + — ^V«i;"^ - 2 j dt j d^^-'^xr^Tr'^^^^Np. (8) 

The symbol stands for the extrinsic curvature of Bt as a surface embedded in S^, 

= (n" — sinhr? m") and is the versor normal to Bt in Sj. The subscript a in Trf^ 

cosh T] 

is a reminder that it has to be considered a tensor density with respect to the measure ^/a. 



3 The field equations 

3.1 The conjugate momenta 

Variation of the action with respect to A^* gives 



where V is the covariant derivative with respect to the induced metric Qij. In a York 
K = const, gauge we have 

Tj 

2V^V,^ = -y^S'ix - qn)Pni (10) 

where vr^ is the traceless part of ti'^K When K vanishes, tt*-' becomes traceless since in 2 + 1 
dimensions 

Ti'^ ^^{K'^ -g'^K) (11) 

and thus vr^^ has only two independent components. Using the complex coordinates z = x+iy 
and z = X — iy, eq.(|TO|) reduces to 

n 

5.vr\- = -i5^P„,-52(^-z„) (13) 

n 

with 7r^2 = vr^^ — i ir^y and J 6'^{z)dxdy = 1. They are solved by 

n\ = -—y^-^ (14) 

n " 



Zir ^-^ z — Zn 
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We shall work in the cm. frame i.e. with V P„ = 0. With such a restriction tt^, and vr^- 
decrease at infinity at least like l/l^;]^. In principle one could add to the solution (0) (and 
(p^D) and arbitrary analytic (antianalytic) function. Our choice to maintain the behavior 
tt\ = 0{l/\z\'^) will give rise, as discussed in sect. 3. 3 to a well defined asymptotic behavior 
of the conformal factor e^'^ ~ {zz)~'^ and ~ ln{zz). Thus the choice of solutions (|n|,|15|) 
is a way to impose the form of the asymptotic metric. 

As we already mentioned in the introduction and will be shown below, one can apply the 
instantaneous York gauge K = only to the topology of the plane (open universes) or the 
topology of the sphere. The sphere can be described by the metric 

where a is a conformal factor regular at infinity. The equations for the tt*^^ are still provided by 
eq.(P,|l^) and their solution by eq. ([l4| , [T5D . On the other hand for the sphere the point z = oo 
is a regular point. Imposition of the regularity of the tensor g^^^'^7i% at infinity, obtained 
through the transformation z' = 1/z imposes for vr^^ the behavior vr^^ ^ z"" with n > 4. 
This is refiected on the following sum rules on the momenta J2n^riz = 0, J^n^nzZn = 0, 
Yln-^^-z^n — 0- Thus except for the static case in which all momenta vanish we need at 
least four particles on the sphere. In the following we shall be mainly interested in the plane 
topology. 

The tracelessness of the momenta vr*-' in the K = gauge entails some conservation 
laws which relate combinations of particle positions and momenta with no reference to the 
field variables. We shall see in sect. 7 that such relations can be derived directly form the 
equations of motions. 

We want to examine here the expression of the conservation of angular momentum in the 
K = gauge in the context of the hamiltonian formalism. Under a rotation 

Qijix^) gij{x^ - ex^ei^) - et/ g.yj{x) - etf gij\x) (17) 
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and similar transformations on the remaining field variables, the volume hamiltonian ex- 
tended to a disk of radius r and the boundary hamiltonian computed on a circle of radius 
r are left invariant. The term ^„-Pni5'^ is also left invariant while the term n'^^gij is left 
identically zero. As a consequence the quantity 

(18) 

n 

is a constant of motion. This is the expression of the conservation of angular momentum in 
the K = gauge. 



3.2 The conformal factor 



The variation of the action under gives the hamiltonian constraint [jTB|, [171 



H = — 



vr fevr „ - (vr 



c \2 



y^R + J2 ^^(^ - ^n) \/ml + g'^Pn, = 0. (19) 

n 

which in the instantaneous York gauge, i.e. K = becomes the non linear equation 

2Aa = -7r>V-'" - Yl ^'^^ " ^n) + 4P„,P„, e'^-. (20) 



Integrating eq.(|TP|) over all space and using the Gauss-Bonnet theorem we see that, due to 
the positivity of vr'^^vr^^ the K = gauge is applicable only to closed surfaces of genus or 
to open universes M. We shall be concerned in the following with open universes. 



The term vr 



a 



27r^^7r^g is the product of an analytic tt^^{z), and an antianalytic 



function and as such the laplacian of its logarithm is zero except on the singularities. The 
structure of vr^, in the cm. frame where J2n^n = 0, is 



IT 



(21) 



with pj\j-_2{z) a polynomial of degree Af — 2, being Af the number of particles. It will possess 
M —2 zeros in the complex plane and this is the origin in the present approach of the so called 
apparent singularities. Thus in the ADM approach the apparent singularities are given by 
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the zeros of the extrinsic curvature tensor Kij because we have the general relation eq.(|TT]) 
and in our gauge K = 0. Due to the Gauss- Codazzi relations they are also the points where 
the intrinsic curvature of the two dimensional time slice vanish. We shall denote by z^^ the 
location of the apparent singularities. We shall now go over to the function 

2d = 2a- ln{2n\TT%) (22) 

and thus we have 

2Aa = -e-^^ -J2s\^-^n)i V< + 4P„,P„,e-2- 4^'^' " ^a) (23) 

n A 

which we want now to discuss near the particle singularities. 

Near the n particle, which to simphfy the notation we suppose to be placed at z = 0, 
eq.(^) reduces to 

2Aa = - 47r(a„, - l)S\z) (24) 

with 

47ra„ = + 4 CnPn.Pnz (25) 

and Cn = 6"^'^*-°''. Cn cannot be infinity otherwise a is not a solution of eg . (|2^) . Eq.(^) can 
be solved exactly |19| by performing the transformation w = [z/AY" on the solution of the 
Liouville equation 2A(x = — e"^*^ representing the Poincare pseudosphere, to obtain 

e-"' = ^ ./~7 .2 (26) 



lal {zz/A^Y"-^ 
A^(l - {zz/A^Y"f 



with arbitrary A. In particular, for \z\ ^ we have 

,2 / s a„-l 



6az (zz 



Combining such a behavior with the factor vr^^vr''^ and keeping in mind that m„/47r > we 
have that the factor e~^'^ vanishes on the sources, which implies 47ra„ = m„. This simplifies 
the source term in eq. (p3D . We have 
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The analogous calculation on the apparent singularites gives a e ^'^ regular at those points 
as expected from eq. (|20|) and the equation for 2cr takes the simpler form 

2Aa = - 47r ^ 5^{z - - 1) - 47r ^ 5^{z - za). (29) 

n A 

3.3 The lapse function 



The variation of the action with respect to Qij [IB, [T^ gives 



p p ai bj 

+V^(7r'W'") - vr^^V^iV^ - vr^^V^iV^- + iV V ^^(x - g„) ^ "^f ^ == (30) 

where we took into account that in two dimensions R^^ — = 0, and the variation of 

the action with respect to vr*-' gives 

2N 

g^J = —{TTij - QiXa) + VjN., + \/^N,. (31) 

Using the time evolution given by the above two equations, the relation K = const = 
becomes 

AN = n%7r\e-^''N + iVe-^- V ^^J^s^z - z„) (32) 

^ 2m„ 

n " 

where we took into account that on the particles e~^'^ vanishes. This is a linear partial differ- 
ential equation for N. We are interested in the behavior of N near the particle singularity. 
Near the particle n the equation takes the form 

^ PnzPn2 2aj^ ^3l£3ll e''^ N 6\z) . (33) 

27r^zz 2m„ ^ ^ ^ ^ 

A solution to eq. ( P3| ) must be such c = lim^^o Ne~'^"' is finite. Moreover it is not difficult to 
prove that c = 0. In fact by integrating twice eq. (|33|) in \z\ we find 

N ^ In^ + ci In 1^1 + Co, (34) 
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which substituted in the definition of c gives c = 0. Moreover as for c = the 6 contribution 
in eq.(^) is absent we have also Ci = 0. The conclusion is that is finite at the particles 
and the S term in eq.(^) is absent due to the vanishing of A^e~^°" on the particles and 
satisfies the linear homogeneous equation 

AN = e-^^N. (35) 

Now we analyze the equation 

AN = ^Jl^e-^-N (36) 

for small z. Substituting the behavior of e^^"^ we have 

A^-f (|)"""iV. ,37) 

Going over to the variable r = \z\/K and integrating twice we have for the s-wave around 
the singularity 

l^r^N = %^^l{rY"-'N (38) 
r dr dr 

and we find 

N = const. [1 + 2{zz/A'^Y- + 0{z) + 0{z)]. (39) 
We examine now the large distance behavior of the conformal factor. Integrating 

A(2a) = -n%7r\e-^'^ -J^S'iz- z^)m^, (40) 

n 

we obtain 

Urn 2 (p Va ■ n dl = -M = -Ann, (41) 



where 



M = ^m„ + j 7r>^„e-2" d^z . 



(42) 
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which gives for a the asymptotic behavior 

M 

2(7 — ln(;z^) + In s. (43) 

This result holds if the integral appearing in eq.(^2|) converges. As ^„-Pn = we have 
TT"'i,TT^a ~ integral converges provided — 1 + /i < 0, which means that the opening 

of the cone at infinity cannot be negative. (We notice that with a choice Pnz 7^ 
the integral in eq.(|^) would be divergent). Actually one has infinite solutions to eg. (^0]) 
depending on the total energy M of the system which has to be specified in solving it. We 
shall see at the end of this section that M cannot depend on time. 
Similarly applying Gauss theorem to 

AA^ = e-^^N = TT%'iT\e-^'' N (44) 

we obtain 

71 

N Inizz) (45) 

An 

with 

n = I d^zn^n^e-^'^N. (46) 



The behavior eq.(^) is consistent with the convergence of the integral eq.(|4B|) provided 
M/An = < 1 which is the condition we already met. On the other hand it is easily seen 
that a behavior A^ ~ \z\^ which makes the integral (^) divergent at infinity is inconsistent 
with the differential equation (^^ . 

We want now to connect the solution 2a of eq.(PD|) with the solutions A^ of eq. (P^D . We 



already noticed that 2cr depends on the free parameter M, which is the total energy of the 
system, and such a parameter does not appear in eq. (pOf) . Then it is immediately seen that 



N = ?t^ (47) 
is a solution to eq . (|35D . The equation for A^ being homogeneous, its normalization is con- 



ventional and could also be taken time dependent. With our choice eq.(H^), the behavior of 
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N at infinity is, from eg. (^3]) ~ ^ln(z^). On tlie otlier liand it is very simple to prove 
tlie following uniqueness theorem: If behaves at infinity like ~ Cq + Cir~", a > and 
on the particle singularities as in eg. (P5|) , then = 0. 
In fact multiplying by A^ eg . (|35[) 



/ d^z V{NVN) = [ d^zVN-VN+ [ d^z e-^^N^. (48) 

J A J A J A 

But for r — > oo 

A^ VA^ -ndl^ 2n{-a)r^'' (49) 



and on the particle singularity for r — i>0 

/ N VN ■ndl^2nr[co + Ci{r^Y-][2cirfin{r^f^'^]^0 (50) 

JCr 

which proves our assertion. 

3.4 The shift functions 

The traceless part of eg. (|3l|) gives 

d,N' = -7r% e-^'^N (51) 

d,N' = -Ti\ e-'^AT. (52) 
Use eg.(|4^ on the r.h.s. of the first. Then we have 

2TT\{z)d-,N' = -4d,d,N (53) 

whose general solution is 

= -^--d^N + g{z) (54) 



13 



where g{z) is an analytic function of z. A similar solution holds for A^^. The analytic function 
g{z) must be chosen as to kill the poles generated by the zeros of tt^^{z), which occur only in 
presence of three or more particles, and if we are interested in describing a reference system 
which does not rotate at infinity g{z) has to be chosen as to give A^^ < \z\ at infinity. We 
shall see in sect. 7 that, due to Xln-^nz — O5 such a function can diverge at infinity at most 
linearly. As we shall see g{z) encodes the time evolution of all the quantities of the problem. 
The trace part of eq.(|3TD gives 

2& = N'd,{2a) + d,N' + C.C. (55) 

By substituting the asymptotic behavior for 2a eq.(|^) into eqs.(|^,^) we see that M = 
as the logarithmic term M\n{zz)/A7c on the l.h.s. is not matched on the r.h.s.. 



4 The particle equations of motion 

In the present section we shall derive the general expression for the time derivative of the 
particle position z„ and of the particle conjugate momenta P„. 
Variation of the action with respect to Pm gives 

z^ = -N\zn,t) = -g{zn). (56) 

The equation for P„ 

d (Pr.aN'^ - N^P^aPubg^^' + ml) 
Pm = Q- H^n) (57) 

requires a little attention because in general the metric is divergent at the particle position. 
Let us consider the expression of Pn 

P„.^.m„^MM&±£^^ (58) 



"n J 

^Jm-g,,{zl + N%zi + N^) 



At the point z = Zn we have + A^-' = but it does not mean that P„j = as gij = gijC^" 
is divergent at that point. We have to compute the behavior of N^{z) — N^{zn) for z Zn- 
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From eq. (|5^ ) we have 

iv.(,) _ « H!!^ ( ^.-z^-z„) y ^ _^ ^^^^ j^gj 

and 

e^<^^^^^ ^ T7— T +0{z-z^) + 0{-z--z^). (60) 

The square root in the denominator goes over to N{zn) and thus in the hmit 2; — > z„ we have 

m..^^sMi±£^ P„,. (61) 
-9„(ii + N')(i'n + Ni) 

With regard to the convergence of the numerator we notice that it is of the type 

{{Z - Zn){z - Zn))-^-\{{z - Zn){z - Z^))^" + 0{z - Z^) + 0{z - (62) 

which converges absolutely to 1 only for /i„ < 1/2 i.e. for not too heavy particle. On the 
other hand if we average the l.h.s. of eq.(|6TD on a circle around Zn and then take the limit 
for zero radius it always converges to Pnz- We find a similar problem in writing the equation 
for Pnz- Taking into account the behavior eqs.( |5P| , |55D of the functions and A^^ for z ^ Zn, 
eq.(0) becomes 



P -P ^-m ^ (63) 

OZ OZ 

and a similar equation for Pnz- Eq . (|63D has to be understood as taking the average over a 
circle around z„ and then taking the limit for zero radius. The procedure is similar to the 
one adopted in |T^. The situation is analogous to the one which occurs in electrostatics in 



presence of point charges. The electric field is infinite or better not defined on the particle, 
but the average of the electric field on a small sphere surrounding the particle gives the 
effective field which acts on the particle as in such a process the field generated by the 
particle itself averages to zero and what survives is the external field. We shall see several 
application of such equations in the following. 
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5 The two body problem 
5.1 The conformal factor 

In the two particle case the equation for a has only two sources. 

A(2a) = -e"'^ + 47r(l - fii)6\z - z^) + 47r(l - fi2)6\z - z-,) (64) 

and we have no apparent singularities as t^^-Xz^ has no zeros in the complex plane. The 
procedure for solving such equation is well known ||19[ 



2^ 8f'(z)f'(z) , , 

(1 - ]\z)]Kz)? 

where the function is given by the ratio of two independent solutions of the second 
order fuchsian differential equation with two regular singularities in z\ and Z2 and one at 
infinity pi 



In fact neglecting for the moment the sources one has simply to check that, being 
ln(/'(z)/'(z)) an harmonic function, 



m,d^\^iX-j]f = -^-l—. (67) 



On the singularities the difference of the indices given by 

y.^iz-Zn)"", y2^{z-Znf- (68) 

has to be such that on the singularity the term ln(/'(z)/'(z)) gives 

49^<9g(a„ -/?„-!) ln((2; - Zn){z - Zn)) = -47r(l - yU„)5^(2; - Zn) (69) 
i.e. an — Pn = fJ'u- The differential equation for y is provided by 

2/"(C) + Q(CMC) = o (70) 
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with ( = {z — Zi)/{z2 — zi) and 

^^^^"4 1 +(C-i)^ C(i-C) J- ^ ^ 

In fact with two singularities at the finite and one at infinity the hnear residues are fixed 
by the Fuchs relations in terms of the /ii, /i2 and /ioo (see also eq.( p.27| )). To understand the 
behavior at infinity let us rewrite 

(2/2^/2 - kk yiyiY 

where W12 = k{y[y2 — ^1^/2) is the constant wronskian. From eq.([T^) we see that at infinity 

_2<T^-(l+/Xoo)ln(CC) (73) 

or 

_2a = -2a + Htt^tt',) ^ ^ In(CC) - 2 In(CC) (74) 
i.e. /ioo = 1 - M/4n = 1 - /i. 

A further requirement on the mapping function / is to give a well defined, i.e. monodromic 
a when one goes around the singularities 0, 1 and consequently 00 in the variable (. From 
the general theory of fuchsian differential equations we know that the vector {kyi,y2) goes 
over to a SL{2,C) transformed vector {ayi + by2,cyi + dy2) due to the constancy of the 
wronskian. From eq. (^) we see that in order that e'"^" be well defined we need 



a h 
. c d, 



eSU{l,l). (75) 



The simplest way to realize this constraint is to choose / = kyi/y2 with yi and y2 behaving 
as in eq . (|68|) ; then / satisfies the monodromicity condition around the point C = fo^' 
k. The point now is to determine k such as the a is monodromic also around the point 1 
and as a consequence at 00. For completeness we give in the Appendix the expression of yi, 
j/2 and the derivation of k. 
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The so determined value of k gives through eq.(|72D and eq.(^) the conformal factor e^°". 
Then we obtain from 

whose behavior for large \z\ we already know to be ~ j- \n{zz). 



5.2 Determination of the trajectory 

It is interesting that in order to write down and solve the equations of motion of the particles 
we do not need to know the explicit form of 2a and A^. In fact from eq. (^6]) we see that 
what we need is in zi and Z2. For two particles let 

2P = Pi- P2; P = Pi (77) 



from which 



and thus 



and 



271 {z - Zi){z - Z2) 



a h ^ ^PzPz {Zi - Z2){Zi - Z2) 



ATI (Z — ZI)(Z — Z2) / N / N 

= — -d.N + g{z). (80) 

-< 2 ^2 — Zi 



The first term at infinity behaves like 



z 



1) 



and in order to have a non rotating reference frame at infinity we must choose g{z) such as 
to remove this divergent behavior at infinity i.e. 
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As 7r^2 has poles in zi, we have 



and then defining I = z\ — z^^q have 



1 = ^. (84) 



We now turn to the equation for P^. We show in general i.e. for any number of particles 
that is given in terms of the residue of the simple poles at the point z„ of the function 
Qi^z^ which appears in the fuchsian differential equation. 
Let the singularity at Zn be 

Q(^) = ^^ + ^ + Qn(^) (85) 

where for the time being z denotes the difference z — Zn- The expansion of the singular 
solution 2/2 around the singularity is 

y2 = z^{\^h^z\0{^z'')) (86) 

with 



2(1 -^„) 

Then substituting into eq.(|63D we have in the limit 2; — > ^„ 



dN°- dN ddr, 

Pn. = Pna-T- ' = + ^n.^?'(^n). (88) 

oz oz oM 

As we already mentioned, such formulas hold in the general A/'-particle system. In our 
specific two particle case we have from eq. ([7T|) 
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where working in the z variable the /3's get divided by 2:2 — with respect to the /3's appearing 



in eq. QTTl) . Substituting into eq. (|88D we have 



47r(Z2 - Zi) Z2 - Zi 

Thus we have to solve the elementary system 

/ = ^; A = -f (91) 

Taking the time derivative of the first equation and introducing the variable A = — we 

at 

reach 

^ = (/X - l)dt (92) 

which solved gives 

/ = const [(1 to) -^^/2]^ (93) 

which agrees with the solution given in 0. 

It is interesting that the same solution can be obtained using simply the first of eq.(pT|) 
i.e. / = I/-P2 and the general conservation laws we shall derive in sect. 7. For two particles 
the conserved angular momentum is given by 

L=Y, xle/Pnj = tm - IPz)- (94) 

n=l,2 

On the other hand eq. (|131|) gives for two particles 

IP, + IP, = 2(1 - fi)t + const (95) 

which gives 



IP, = (1 - /i)t - + const. (96) 



Using now I = l/P, we reach immediately the solution eq. ( 
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6 The asymptotic metric 



In the previous section we gave the exact metric for the two particle case. In case of M 
particles again the solution of the Liouville equation is given by eq.(|72|) where y2 are 
solution of a fuchsian differential equation with singularity at infinity of index /Xoo = 1 — 
M/47r = 1 — /i. Then the asymptotic behavior of a for large z is given by 

2a = 2 Hy^m - k^ym) - In |m;i2|' ~ (2 - /i) Hzz) - k^zzY'^ + ... (97) 

Thus 

e^'" ^ const {zzY^". (98) 
It follows that the asymptotic behavior of is 

^ = ^fe^ ^ i ^^^"^^ + ^ ln(z^)(zz)^~^ + ... (99) 

and then 

= const z{l + k'^ifi - 1) \n{zz){zzY-^ + ....) + g{z). (100) 

If g{z) is chosen as to remove the behavior proportional to z (non- rotating frame at infinity) 
the behavior of becomes 

A^^ ^ const z\n{zz){zzY-\ (101) 

, p.01| ) provide the general asymptotic behavior of the metric. 
In this formalism the appearance of a Gott pair is signalled by an oscillating behavior of 
f{z) as \z\ oo (the total mass is now an imaginary number) 0. This pathological behavior 
is not surprising because both the formalism and the gauge used here assume the existence 
of a global time. 



Formulas 



7 The A/^-particle case 

We shall consider in this section the A/'-particle case with particular attention to A/" = 3. 
Again the key role is played by the conformal factor e^°" or equivalently e^°" which is given 



21 



by eq.(^) where yi, 1/2 are solution of a fuchsian differential equation. It is interesting that 
can be rewritten in the following form 

^^0Gm^2jaj)^^^^^^ ,102) 



where 



dM 

In fact 



Gif) = Mj- (103) 



dM f'dM 1 - // 

and it is easily checked that 

I df d [ df 



d(-2a) 1 df 2fG(f) , 
^ - + + c.c. (104) 



f dM df \dM 
With regard to the conjugate momenta we have 



(105) 



where as Yln^ri ~ 0, p^f-.2{z) is a polynomial of degree Af — 2. In flat space the number 
of Lorentz invariants is easily seen to be SA/" — 3 as for n > 3 the 2 + 1-vector P„ is 
determined by the scalar product of P„ with Pi, P2, P3. In curved space we have the same 
number of invariants. As it is well known here the invariants are replaced by the traces 
of the holonomies around the world lines of an arbitrary collection of particles [|l], |21|, p2 |. 



The parallel transport monodromy matrices which belong to 5*0(2, 1) or in the fundamental 
representation to SU{1,1) are defined up to a conjugation. Keeping in mind that every 
SU{1, 1) holonomy has three real degrees of freedom we have 3A/' degrees of freedom to 
which we have to subtract the three degrees of freedom of the SU{1,1) conjugation thus 
reaching as expected the same number of invariants as in the fiat case. 

Given at the time t = the position z„ and the momenta P„ of the particles we know 
from eq.(|n|) also the position of the apparent singularities za- Knowledge of the masses m„ 
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fix also the second order residues of tlie poles in Q{z) at z„ to (1 — fin)/4: while the second 
order residue at the apparent singularities in fixed to —3/4 [§. 

In order to understand the problem better let us consider the case TV = 3. It is useful to 
go over from the equation in projectively canonical form 



y" + Qy = o 



where the general form of Q{z) is p|, ^ 



1 f^n ^ f^n 



to the equivalent form M 



B 



+ 



B 



A{z - zbY 2{z - zb 



y" + py' + qy = 



with 



Q = q 



I 1 
p p 



2 4 

The transformation can be so chosen as to have the following Riemann scheme 

f 1 z^ Za oo ^ 
Poo 

/i2 /is 2 Poo + /ioo/ 

where the Fuchs relation 

/il + /i2 + /is + /^oo + 2poo = 1 



fixes Poo in terms of the other parameters and 

p{z) 



1- fil , 1 - /i2 _^ 1 - P3 



z — 1 z — Z-i z — Za 



(107) 



'10^ 



(109) 



;iio) 



;iii) 



(112) 



(113) 



Z3{Z3-1)H ZA{ZA-I)bj 



z{z-l) z{z-l){z-Z3) z{z-1){z-za) 



(114) 
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with K = poo(Poo + Poo)- Here we use the unconstrained parameterization of ref.^; the 
connection with the parameters p„ and /3„ can be easily obtained by comparing the residues 
in eq.(|llUD. 



The parameters /ij and poo ai'e free real parameters which give the trace of the holonomies 
around the particles and infinity (total energy), za is given in terms of the and the P„ 
{n = 1,2,3) while H is given by the no- logarithm condition i.e. the requirement that the 
point Za is a regular point for 2a [^. Such H is easily found 

Za{za - 1){ZA - Z3) / 2 //^l , Us -I K \ 
H= rr \hA-{— + 7 + )hA + —, rr > • 115 

^3(2:3 - 1) Za za-1 za- z^ za[za - 1) J 

Thus we have the following free real parameters pi, p2, Ps, Poo, R-e(zA), Ini(-2yi), Re(6A), Ini(6A) 
which are eight real parameters, to be contrasted with the six real invariants we have enumer- 
ated at the beginning of this section. The reason is that the constraint that the holonomies 
described by eq.( |109| ) belong to S'f/(1,1), fixes the value of Re(6A), Ini(6yi) in terms of all 
other parameters. The explicit form of this relation has been investigated in the mathemat- 
ical literature but we are not aware of any explicit form of it. Adding a new particle 
introduces a new particle singularity with exponent p4 i.e. p{z) acquires a new term of the 
type 

(116) 

z — Zi 

but in addition a new apparent singularity is generated say at zb- The indices of this 
new singularity have to be integers, and 2 while the new simple residue hs has to be a 
function of all other parameters, and in particular of zb-, because the SU{1, 1) character 
of the holonomies allow only an increase of 3 in the number of free real parameters i.e. 
P4, Re^;^, Imz^, but as we said above we do not know an explicit solution of this problem. 
Adding other particles does not alter the procedure. 

We come now to the problem of determining the change in time of the position of the 
particles z„ and of the position of the apparent singularities zb- If f{z) = provides a 

monodromic 2cr, N is still given by eq. (^) as the sources are independent of M; the equation 
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for , are solved as before by 

The new fact with respect to the two particle case is that now 'K^zi^) has J\f — 2 zeros in the 
complex plane which are the position of the apparent singularities. We write 

^M^) ^ I ^ Pnz ^ rigl^-^i?) ..-.ox 
2 ATT^Z-Zn V{Z) ' ^ ' 



n 



Taking the derivative of the previous equation with respect to z and then setting 2; = 2;^ we 
obtain a relation which will be useful in the following 

which for a single apparent singularity reduces to 

The general expression of is 

N^= f^^^ , d,N + g{z) (121) 
[[b[z - Zb) 

and thus to remove the poles from we must have 

^dpB 1 V{zb) ^ f . 



where pi{z) is a polynomial, as 



dM^B) = (123) 



In fact the expansion of the singular solution around the apparent singularity za is given by 

y2={z- Za)-H^ + - ^a) + 0{{Z - ZA?)) (124) 
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which substituted in 2a gives the above result. Recalhng the equation of motion of Zn we 
have now 

Thus the time variation of z„ is given in terms of the momenta P„ which determine completely 
the polynomial V{z) and the derivative with respect to M of the residues (3a- 

We come now to the equations for P„. We already saw on general grounds that 

Pnz = 471"^ + Pnz 9 (Zn) = 

= An^-P + P n'(z) (126) 

From the fuchsianity of the differential equation in projectively canonical form we know that 
5^/?„ + = 0; 1 - /iL = 5^(1 - /i' + 2/3„2:„) + J](-3 + 2/3b^b) (127) 

n b n B 

so that due to /Xoo = 1 — M/47r, 

E5M+>^aM=°' 4^(l-4;^) = 2:aM"« + 2:aM^- (128) 

n B n B 

We see that provided pi{z) is a first degree polynomial, p{z) = cq + Ciz, the consistency of 
the equation J^n^n = is assured by the the first equation in (|128|) . Similarly one easily 
checks that the second relation in eg .( [1281 ) provides the law 

n n 

where we have exploited eq.( |119| ) and Pm/ {zb — z^) = 0. Such equation gives both the 
law of conservation of angular momentum 
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which we derived in sect. 3 and the "generahzed conservation law'' 

d sr-^, - f> \ ^ \ 



dt 



(131) 



which appears to be related to a dilatation transformation ||2^ . To complete the scheme we 
must provide the expression of the time derivative of the position of the apparent singularities 
za and the time derivatives of the simple residues at such singularities (3a- The equation for 
za is easy to obtain. One takes the time derivative of 



= 5^ 



Za ~ Zn 



(132) 



which allows to express za in terms of Zn and Pm which are already known. Using eq.( |132| ) 

df3n 



and eq. (|119|) we obtain 




^ [zA + CiZA + Co] + 2^ 



1 



dM Za — Zn 



dl3_ 



B 



1 



DM Za — zb 



n 



C^A\^A ~ ^Cj_^ Ilc^B^^B - Zc) 



(133) 



Using the no-log condition the second term in the above equation can be rewritten as 



-A 



B 



dp. 



dM ^^dMzA- ZB 



(134) 



and thus we reach 



Za = -Co - ciZa + (3, 



d(3. 



_ \- d(3B 1 

dM Ylc^A^ZA - zc) ^ dM 

Za Zb Wci^B^^B Zc) 



V{zb) 



(135) 



In the case of three particles i.e. one apparent singularity, the above equation reduces to 

d^A 



Za = -Co - ciZA + V{za)I3a 



dM' 



(136) 



In order to connect with the mathematical literature we shall choose the meromorphic func- 
tion g{z) such as to keep the position of particle 1 at and of particle 2 at 1. For clarity we 
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shall refer first to the case of three particles i.e. one apparent singularity, and then give the 
generalization to M particles. Thus 



dp. 



z{z-l) 



dM {z - za)za{za - 1) 
Then substituting in eq.(|136D we reach 



1 



dM ^z — za ' za ' za(za — 1 



+ — + 



(137) 



Za = ^(2;a)t77t ( Pa - 



1 

Za 



dM V za - 1 

If we divide the previous expression by is which is simply given by —(7(2:3) we obtain 



:i38) 



dzA 
dz?. 



dH ZA{zA~l)izA- Z3) 



dh. 



Z^{,Z3 - 1) 



Za 



1^2 



Za — ^ Za — Z3 



(139) 



i.e. the first Garnier equation, with hamiltonian eq.( |115| ) where we took into account that 

(140) 



Pa , 1 n - ^1 1 - /i2 , 1 - /is 

-TT = "A + t: [ ^ 7 

2 2 \ Za za-1 Za- z^ 

We come now to the more difficult task of finding the time variation of the simple residues 
Pb- To accomplish this we shall exploit the equation 

d{2a) 



dt 



(141) 



evaluated in the neighborhood of an apparent singularity. The above equation embodies 
part of the information carried by the dynamical equation for gij. 

It is more useful to rewrite equation ( |141| ) in terms of the reduced conformal factor 
2a = 2a — ln(2'7r^j7r^^). In order to accomplish this we need the evolution equation for vr^^ 
and TT^j i.e. for the traceless part of vr"^. These are given, outside the particle singularities 

by 



(142) 



Using eqs. (|55|,p!42D we reach 



2d = g{z)d,{2d) - g'{z) + c.c. 



(143) 
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We now recall the fundamental relation of the theory of fuchsian differential equation |T9[, 
which can be easily proved starting from eq. (^) which holds in the general A/'-particle case 

Q{z) = -\[dl{2a) + \{d,{2~a)f] (144) 

where the r.h.s. plays the role of the analytic component of the energy momentum tensor 
in two dimensional conformal Liouville theories. Taking the time derivative of eq. ( |144D we 
have 

Q{z) = \g"\z) + 2g\z)Q{z) + g{z)Q\z). (145) 

It is interesting that this equation contains the whole system of the Garnier equations for the 
time evolution of the apparent singularities. If we Laurent-expand g{z) about the apparent 
singularity za 

g{z) = + g, + g,[z- z^) + ^{Z- Za)' + 0{{Z - ZA?) (146) 

z — za 2 



and 



Q{^) = -TT^-^ + ^TT^^ + Qa{za) + Q'a{za){z - za) + 0{{z - za)') 
4:[z — Za) Z[z — Za) 



+ wr^^ - ^ + Q'a{^a){z - za) + 0{{z - za)') (147) 



4{z - za)^ 2{z - za) 4 

The fourth order poles in eq. (|145|) cancel identically without producing any information; the 
matching of the residue of the third order pole gives 

Za = 9-iPa - 9o (148) 
while the matching of the first order pole gives 

Pa = -2g^iQ'A{zA) + gil3A - \g2. (149) 
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Eq.( |148| ) reproduces the results obtained in eq.( [135| ), as can be seen by computing 
qq form eq.(p.22|). As an example we can apply the above equations to the three 



case. If we keep particle 1 at an particle 2 at 1, the g{z) function takes the form 



DM 



+ 



1 



Z- Za ZAiZA - 1) za-1 



1 

Za 



Then substituting into eq. 



Za = V{za 



we have 

df^A 



dM 



/9/ 



Za 



1 

Za 



zz 



dH 



and 



^A = V{ZA 



dl3A 
^ dM 



'^Q'a{za) 



-zs 



dH 

dzA ' 



'za - 1 

These are the Garnier equations with hamiltonian eq.( |115D once we keep in mind 
cording to the general equation eq.( |125| ) 

^3- ^'^>dMzAizA~l)izA-Zsy 

It is important to remark that the matching of the residues of eq. ( |145D on the 
singularities reproduce the equation of motion z = —g{zn). 



g-i and 
particle 

(150) 

(151) 

(152) 
that ac- 

(153) 
particle 



8 Conclusions 

It is useful to summarize here the basic features of our derivations which are rather straight- 
forward. In the instantaneous York gauge K = the momenta vr"^ conjugate to the space 
metric have only two independent components vr^^ (^r^^) which are meromorphic (antimero- 
morphic) functions of z whose residues are the particle momenta P„. Knowledge of such 
momenta 71% allows us to write a Liouville equation for the reduced conformal factor a, in 
which the sources are the particle singularities i.e. the poles of tt^^ and the apparent singu- 
larities, i.e. the zeros of vr^^. The function is given by the derivative of a with respect to 
the total energy while is given in terms of by 

= -^d,N + g{z), (154) 
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where g{z) is a meromorphic function which cancels the polar singularities of the first term 
and grows at infinity not faster than z. The linear term in g{z) is fixed if we want to deal 
with a reference frame which does not rotate at infinity. 

The equation of motion for the particle position are i„ = —N^{zn) = —g{zn) while those 
for the particle momenta 

Pnz = PnadzN'' - m^d^N. (155) 

On the other hand the change in time of the position of the apparent singularities and their 
residues are given by the ADM equation for a computed on the apparent singularities, thus 
providing the Garnier hamiltonian system. 

In the solution for and consequently in g{z) and P„ there appears the derivative with 
respect to the total energy M of the coefficients of the fuchsian differential equation which 
underlies the determination of a. I.e. the change in the /3's has to be such as to remain in 
the SU{1, 1) class of differential equation. It has a very simple solution in the case of two 
particles while in general it is obviously related to the Riemann-Hilbert problem even if it is 
less specific. 
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The two solutions of eq.(|TOD which near the origin are a power multiplied by an analytic 
function are given by p5 



2/.(C) = C^(C-l)^«.(C) (156) 

where 

ui = F(a,6;c;C); U2 = C}^'F{a - c + l,b - c + 1-2 - c-Q (157) 

and 

1 1 

a = -(1 + /ii - /i2 + /Uoo); & = 2(1 + /"i "~ - /ioo); C=l + /ii. (158) 
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Under a circuit around point 1 they mix as Ui BnUi + Bi2U2 and thus the demand that 
{kyi,y2) transforms under a representation of SU{1, 1) imposes 

\k'\ = ll^l (159) 



B 
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We have 25 



D o^^ J-^ic-a-b) r(c)r(c — 1) nfin^ 

r(c — a)r(c — o)r(o)r(a) 



r(i-a)r(i-6)r(i + a-c)r(i + 6-c) 



(161) 



which gives 

\k^\ = |A(a)A(6)A(l - c)A(2 - c)A(c - a)A(c - b)\ (162) 
where as usual A(x) = r(x)/r(l — x). 

(163) 
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